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Chiral perturbation theory, dispersion relations and final state 
interactions in K ^ tttt 

Gilberto Colangelo^ 

'^Institut fiir Theoretische Physik der Universitat Zurich 
Winterthurerstr. 190, 8057 Ziirich, Switzerland 

We discuss the recent literature on the treatment of final state interactions in K ^ nn. Various approaches are 
compared and particular emphasis is given to the possibility of combining dispersive methods with lattice input. 
Recent results on the dependence of various quantities on the quark masses at order in the chiral expansion 
are presented and the relevance for the lattice calculations is discussed. 



1. INTRODUCTION 

Ideally, lattice QCD calculations should pro- 
vide numbers that can be directly compared to 
the experimental ones. In practice, one often 
needs extra theoretical input and/or treatment 
before being able to compare to the phenomenol- 
ogy. One important example of this is the K 
TTTT amplitude: until the recent work by Lellouch 
and Liischer |l[ it was not even known how one 
could directly calculate this amplitude on the lat- 
tice, and the problems of principle related to a 
naive approach had been identified by Maiani and 
Testa in a famous no-go theorem [^j. The stan- 
dard method to work around this no-go theorem 
consists in calculating on the lattice an unphysical 
matrix element, K n, and to use an approx- 
imate relation (based on tree-level chiral pertur- 
bation theory (CHPT)) to obtain the K tttt 
amplitude ||]. The latter step is known to repre- 
sent a rather crude approximation, and sizeable 
higher order CHPT corrections are expected. Al- 
ternative methods to obtain the K ^ tttt ampli- 
tude were discussed in Despite its theoreti- 
cal cleanliness, the use of the Lellouch-Liischer 
method represents a formidable challenge, and 
is still out of reach for today's state-of-the art 
numerical simulations. Indeed all lattice collab- 
orations 1^ that have presented results for the 
K — > TTTT amplitude at this conference have used 
either the method of Bernard et al. ||^ or one of 
those proposed in [Q. 

In view of this, it is worth to ask whether one 



could devise a better way to obtain the K — > tttt 
amplitude from quantities that are calculable to- 
day on the lattice. For example, since today's 
calculations are all done in the quenched approx- 
imation, one may try to find a way to obtain 
information about the unquenched results from 
(partially) quenched calculations. Work in this 
direction has already been done, see e.g. Ref. [^, 
but I will not touch this issue here. Another very 
important issue is that of final state interactions 
(FSI): if one calculates the if — > tt matrix ele- 
ment and obtains the K — s- tttt amplitude by using 
tree level CHPT, one is neglecting FSI. Although 
these unitarity effects are formally of higher or- 
der in the chiral counting, they are known to be 
particularly large when two pions interact in the 
S wave in the isospin zero channel, and might 
be crucial in the final numerical result. The role 
of FSI has been discussed at length in the recent 
literature, particularly in connection to the calcu- 
lations of e'/e in the Standard Model The 
emphasis of the recent literature is on the dis- 
persive treatment of the FSI, but all the groups 
that calculated e'/e making explicit use of the 
chiral expansion beyond leading order, did also 
find large corrections due to FSI ||l^ , |l3|| . Ear- 
lier treatments of if — > tttt, and in particular of 
the CP conserving part of the amplitude also dis- 
cussed the importance of FSI, and its role for the 
Ai = 1/2 rule |l|,|l|. 

In this talk I will review the recent and the 
less recent literature on the treatment of FSI in 
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K — > TTTT. I will do so by comparing two different 
approaches: one which explicitly relies on CHPT, 
and one that uses only dispersion relations. A 
third possibility, which is also discussed, consists 
in a dispersion-relation inspired improvement of 
a chiral calculation. The focus of this discussion 
will be on the methods, rather than on the nu- 
merical results for the K nn amplitude, which 
I will not touch at all. A second aspect that I 
will discuss in this talk is the dependence of some 
quantities on the quark masses, as it is evaluated 
in CHPT at next-to-next-to-leading order. I will 
present the analyses which are currently available 
for masses, decay constants and scatter- 
ing amplitudes. The latter result is directly con- 
nected to the problem of FSI: it has been shown 
that large FSI are also reflected into a very strong 
dependence on the quark masses [ p^ . 

2. THE SCALAR FORM FACTOR: 
CHPT vs. DISPERSIVE METHODS 

In order to compare the chiral expansion and 
the dispersive approach, it is easier to consider a 
simpler quantity than the K — > tttt amplitude - 
if we want two pions in the S wave and isospin 
zero in the final state, the simplest quantity we 
can think of is the scalar form factor of the pion: 

T{s) =JV{Tr{pi)7T{p2)\uu + dd\0) , (1) 

where s — {pi + P2)'^ and A/" ensures r(0) = 1. 
In CHPT this has been calculated (and studied 
numerically) to two loops and has also been 
analyzed with dispersive methods , so that we 
can immediately compare the two results. We 
first consider the CHPT calculation, which reads, 
to one loop 

r(.) = l+^(4-l)+^(2.-Mj)J»+0(,s2),(2) 

where N = Wtt^F^, and J(s) is a loop function 
||20|| . In order to compare to the dispersive treat- 
ment, it is useful to rewrite this expression in a 
different way, that makes the tttt phase shift ap- 
pear explicitly 

T{s) = 1 + CS + - / ds' } \ +0(.^) , (3) 



where c is a constant related to €4 ~ InM^, and 
where 5^"^^ is the leading order tttt phase shift in 
the S wave and / = channel 

5{s) = 5^^\s) + 5^^\s) + . . . 

- ^(2.-M,^) , (4) 

where cr(s) = \/l — 4M^/s. In order to have a 
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Figure 1. Scalar form factor: comparison of the 
CHPT calculation to leading, next-to-leading and 
next-to-next-to-leading order and the dispersive 
one (see legenda). The solid curve is the "true 
form factor" , see text. 



numerical representation in CHPT to one loop, 
one needs to fix one constant, i^. Once this is 
done, e.g. by fixing the derivative of i^(s) at s = 
at its physical value (see below), the form factor 
is as shown by the dashed line in Fig. 0. The 
same kind of analysis can be repeated to two loops 
- there, however, one has to fix more constants 
|l8|| . In Fig. it is clearly seen that while the 
chiral expansion behaves nicely up to around s = 
(0.4 GeV)^, at higher energies it does not work 
so well. In particular, the behaviour of the series 
at s = Af|. looks as follows: 

1 (LO) ^ 1.62 (NLO) ^ 1.67 (NNLO) , (5) 
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whereas the "true value" is 1.42. At sUghtly 
higher values, the two loop form factor does bend 
down, as the true form factor does - however, at 
a quantitative level it does not give a satisfactory 
representation. It is not difficult to understand 
why two loops are not yet enough to reproduce 
correctly the form factor at these energies. The 
real part of the form factor is proportional to the 
cosine of the phase 



Rer(s) - cos(5(s) - 1 - 



(6) 



and since the phase is a quantity of order in 
the chiral expansion, this negative contribution to 
the real part starts only at NNLO. Numerically, 
at NNLO the cosine becomes different from one, 
but it is still larger than the real value by about 
20 % - which roughly explains the difference seen 
in Fig. |. 

The dispersive treatment for the form factor is 
quite simple: Omnes tells us that if one knows 
the phase cj) of the form factor (on the positive real 
axis, from s = AM"^ up to infinity), one knows 
the form factor in the whole complex plane. Its 
explicit expression reads 



r(s) = exp - / ds -— — 



(7) 



In principle one should worry about the presence 
of zeros, and consider a corresponding polynomial 
factor, but for the sake of simplicity we will dis- 
regard this possibility. We also assume that the 
integral in Eq. converges. 

Note that what enters in the Omnes function is 
the phase of the form factor - according to Wat- 
son's theorem, below the inelastic threshold this 
is equal to the tttt phase shift in the S wave and 
isospin zero. Above that, however, it cannot be 
easily related to observable quantities. Note also 
that although the inelastic threshold is formally 
at s = 16M^, the 4 tt channel is very weak, and 
can as well be neglected. Around 1 GeV the much 
stronger 2 K channel opens up, and one cannot 
neglect inelastic effects anymore. 

In fact what we have indicated until now as 
"the true form factor" is nothing but the result of 
a dispersive treatment which explicitly takes into 
account the contribution of the 2 K channel, in 



a fully consistent manner (coupled channel anal- 
ysis), and using the available data [|l^. In that 
case, one cannot express the form factor as sim- 
ply as in Eq. (0). However, even neglecting the 
inelastic channels, and using the tttt phase shift, 
e.g. as given in [Q, all the way up to a cut-off, 
one gets a decent description of the form factor, 
as shown in Fig. |l|, where the cut off Omnes func- 
tion 



il(s) — exp ■ 



JiMl S'(s'-S) 



(8) 



for A = 1.4 GeV is plotted. The precise value of 
the cut-off does matter for the final result, but 
not too much, as long as one is not taking it too 
low. 

It is interesting to compare at the algebraic 
level CHPT to one loop and the Omnes represen- 
tation: in order to do so we have to expand the 
latter in a chiral series - this is done by expanding 
the phase of the form factor, and correspondingly, 
the exponential. If we do so, we end up exactly 
with Eq. (||). Note that if we simply insert the 
leading order chiral phase, Eq. (^), in the disper- 
sive integral (^ we get a divergent result, and 
therefore are forced to subtract it at least once 
more, as in Eq. (||). Note also that the the dis- 
persive representation explains also the presence 
of a chiral logarithm in the constant c in Eq. (^ : 
with the leading order phase the once-subtracted 
dispersive integral is infrared divergent and ex- 
plodes logarithmically as the pion mass goes to 
zero. 

In summary, in the case of the scalar isoscalar 
form factor at energies around s = M]^ it is 
clear that the chiral expansion is less economical, 
and gives a poorer description than the dispersive 
treatment, even if the latter is done neglecting in- 
elastic effects. This motivates us to apply a dis- 
persive treatment also to the K — > tttt amplitude. 

3. THE K ^TiTi AMPLITUDE 

It is easy to see an analogy between the K — > 
TTTT amplitude and the scalar form factor discussed 
in the previous section, as far as FSI are con- 
cerned. On the other hand, the physical K 
■KIT amplitude does not depend on any kinemati- 
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cal variable, and the dispersive treatment of the 
scalar form factor cannot be immediately applied 
to this amplitude: first we have to extend the def- 
inition of the K TTTT amplitude to make it de- 
pend on one (or more) complex variables. There 
is no unique way to do so, and we will discuss two 
different options. The chiral calculation, however, 
does not need this extension, and the correspond- 
ing treatment of FSI can be discussed directly for 
the physical amplitude. 

3.1. CHPT treatment of FSI 

In CHPT to one loop the K ^ tttt amplitude 
reads as follows 1l5| 



A 



2M|. 



1 

Ml 



{2Ml - Ml) J {Ml 
Ml 



2N 



In- 



1 



(9) 



where A^"^^ is the leading order amplitude, and 
where we have written explicitly only the once 
subtracted loop function generated by the two- 
pion exchange, and the accompanying chiral log. 
The ellipsis stands for other loop functions and 
counterterms. The problem with this amplitude 
is that in order to use it for numerical calculations 
one has to input values for the weak (and strong) 
counterterms, which are largely unknown, espe- 
cially if one wants to identify the contributions 
of individual operators in the weak Hamiltonian. 
Also, the identification of the contribution of a 
certain loop diagram is ambiguous, not only be- 
cause of the presence of the scale ^, but also be- 
cause the break down into individual loop dia- 
grams depends on the parametrization chosen for 
the fields in the effective lagrangian. Note, e.g., 
that logs coming from other loop diagrams have 
to cancel the term proportional to M|- In in 
(^), so that the full amplitude has a well defined 
ruu — rud — *■ limit. Despite these ambiguities, 
one may find instructive to evaluate the pion loop 
contribution in (|^): for /i ~ Mp this contribution 
amounts to about 40 %, which is indeed quite 
large. On this basis, it has often been stated 
that in CHPT the main one-loop correction to the 
K TTTT amplitude comes from tttt FSI | |l5| , p^ , p^ . 
We stress that, because of the ambiguities men- 
tioned above, this statement cannot be made very 



precise. 

Note the analogy to the case of the scalar form 
factor: indeed Eq. ^ looks completely analogous 
to Eq. (H), for s = M^-. However, while in that 
case we could fix the counterterm and the accom- 
panying log (both contained in £4), with the value 
of the scalar radius of the pion, here we do not 
have a corresponding quantity to use as input. 
Numerical statements about this amplitude de- 
pend in a crucial way on a theoretical estimate of 
the low energy constants involved. 

3.2. Dispersive treatment: K off-shell 

In view of the analogy between Eq. (^) and 
Eq. (|) for s = Ml, it is tempting to make the 
momentum squared of the kaon the variable on 
which to construct the dispersion relation. This 
idea has been proposed by Truong p4| , and re- 
cently revived by Pallante and Pich H]. Although 
the idea is seductively simple, it implies a num- 
ber of serious problems, all related to the fact that 
the off-shell extrapolations of a physical hadronic 
amplitude are arbitrary, in contrast to on-shell 
amplitudes, which are unambiguous and well de- 
fined. A discussion of the effects of this arbitrari- 
ness can be found in Ref . , and the interested 
reader is referred to that paper for further de- 
tails. Both papers 0,01 ignore this arbitrariness, 
and choose an unphysical off-shell extrapolation 
for the kaon. In fact both papers quote for the 
"off-shell amplitude" at tree level in CHPT 



(10) 



where A is a constant. Such an extrapolation 
corresponds to using the kaon field in the effec- 
tive lagrangian for going off-shell. This choice is 
unphysical: one could add a term that vanishes 
at the equations of motion (and hence does not 
affect any physical observable), and correspond- 
ingly change the off-shell amplitude ([lo| ) with 
terms proportional to s — M^. There is no way 
to prefer one choice over the other as long as one 
is using the fields in the effective lagrangian: it 
is well known that these do not have any phys- 
ical significance and are simply integration vari- 
ables in the path integral |2^. The calculation 
of the one-loop corrections to (px|), as done in 
H makes the sickness of this extrapolation even 
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more evident: a standard basis of the 0(p'*) coun- 
terterms, as e.g. given in [Q, does not suffice 
to reabsorb the divergences of the one loop di- 
agrams. Off-shell extrapolations can be unam- 
biguously defined by using quark operators, such 
as the axial current or the pseudoscalar density 
- in that case, the off-shell amplitude to leading 
order is different from (p^, and does not change 
if one adds to the lagrangian terms that vanish 
at the equations of motion Q . One- loop correc- 
tions can also be unambiguously calculated. 

In the approach of Refs. p4|j7| introducing the 
off-shell extrapolation is unnecessary, and, as dis- 
cussed above, misleading. At tree level, their pro- 
cedure amounts to multiplying the on-shell tree 
level amplitude by the Omnes function evaluated 
at s = Ml: ^(2) ^ yl(2)f)(M|-). At one loop 
one can extend this procedure in the following 
manner Isl: 



yl(2) + (4) 



(11) 



where A(2^(M|.) is the 0{p'^) part of the Omnes 
function 



4 A/, 



. s'{s' - s) 



(12) 



which needs to be subtracted from the one-loop 
amplitude to avoid double counting. Note that 
since the 0{p'^) phase shift grows linearly with 
s (^), the dispersive integral (|l2|) diverges when 
A — !■ oo - the final result depends on the choice 
of the cut-off. 

While the effect due to FSI is quite substan- 
tial if compared to the tree level chiral ampli- 
tude, the exponentiation of the FSI part of the 
one loop correction gives a much less dramatic 
increase. If we use a cut-off A = 1 (1.3) GeV, 
we get A(2)(M|-) = 0.3 (0.4), and n^M^) = 1.4, 
(1.65). The choice of the cut-off gives an uncer- 
tainty which is comparable to the increase due 
to the exponentiation. Moreover, to obtain a 
numerical estimate from such a calculation, one 
still needs an estimate for the 0{p^) counterterms 
js), as in pure CHPT. The necessity to refer to 
the chiral expansion and to its many constants 
is, in my opinion, the main disadvantage of this 
method. While it is economical and gives a quick 



estimate of the effect of FSI if one refers to tree 
level CHPT, it becomes much more involved if 
one extends it to one loop. In addition, at one 
loop there remain large uncertainties. 

3.3. Dispersive treatment: momentum- 
carrying Hamiltonian 

The necessity to refer to CHPT disappears if 
one chooses a different extrapolation of the phys- 
ical amplitude, and allows the weak Hamiltonian 
to carry momentum, as has been suggested in 
pX| ]. The method goes as follows: consider the 
amplitude 

i=o{MpiHP2)\HK.^mK{qi)) T+{s,t,u){13) 

with the Mandelstam variables s = (pi+P2)^, t = 
ili-PiY, u = (qi-P2)^, related by s + t + u = 
2M^+Ml+q2, where q2 is the momentum carried 
by the weak Hamiltonian. The physical decay 
amplitude is obtained by setting q!^ — {s — M|-, 
t = u = M^). 

The dispersive treatment of a function of three 
complex variable is rather complicated. However, 
the problem simplifies considerably if one neglects 
the contribution of the imaginary parts of D and 
higher waves. With this approximation, the am- 
plitude can be decomposed into a combination of 
functions of a single variable: 

T+{s,t,u) = Alois) + ll[No{t) + 2Ro{t)] 



Ni{t) 



(14) 



where A = Mj. - M^. 

In this approximation, the dispersion relation 
for the full amplitude is transformed into a set 
of coupled dispersion relations of functions of a 
single variables. These can be solved numerically, 
as shown in [0. An important point concerns 
the overall number of subtraction constants which 
is two for the amplitude T+. In order to solve 
the dispersion relations and calculate the K — > 
TTTT amplitude, one has to provide as input the 
two subtraction constants. One of the two can 
be given in terms of the K t: amplitude - a 
soft-pion theorem relates the amplitude at the 
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Figure 2. The function |T+(s,t,u)| plotted vs. 
E = -^s along the line of constant u = : the 
result of our numerical study for different values 
of X are compared to tree level CHPT. 



soft pion point {s = u = Af^, t = M]^) to the 
K TT amplitude up to terms of order M^: 

-^<|j^.„ + | + 0(m;,, (15, 

where N = A^o(M|.) + 2Ro{M'^). Notice that 
although the process involves a kaon, the rela- 
tion is based on the use of the SU{2) symmetry, 
and therefore suffers from corrections of order AI^ 
only. The existence of this soft-pion theorem is 
a crucial advantage of this approach, because it 
simplifies the calculation of one of the two sub- 
traction constants. The key of the problem is how 
to calculate b. This constant is related, e.g., to 
the derivative in s of the amplitude T+ at the soft 
pion point. In a Ward identity that relates 
this derivative to a K ^ n matrix element was 
derived. The calculation of this matrix element 
would determine b. Alternatively, one could try to 
obtain information on the K nn amplitude at 
another unphysical point (i.e. with meson masses 
at their physical values but with the weak Hamil- 
tonian carrying momentum) and use this input to 
constrain the second subtraction constant. 



In the absence of calculations of 6, the nu- 
merical results were illustrated by fixing & at a 
certain value and then varying it within a fairly 
wide range. To fix the central value leading order 
CHPT was used: 

The size of the correction X is at the moment 
unknown, but nothing protects it from being of 
order A/^. In the numerical analysis the generous 
range X = ±30% has been used. The results are 
shown in Fig. |, where |T+(s, Mj^ + M^-s, M^)\ 
versus s is plotted, comparing the numerical solu- 
tion of the dispersion relations to the CHPT lead- 
ing order formula. The latter is what has been 
used so far whenever a number for the K rnr 
matrix element extracted from the lattice has 
been given. The figure shows that large correc- 
tions with respect to leading-order CHPT are to 
be expected. One source of large corrections is the 
Omnes factor due to tttt rescattering in the final 
state. The other potentially dangerous source is 
represented by X, the next-to-leading order cor- 
rection to the relation ( ]T6| ) between a and b. The 
latter could (depending on the sign) in principle 
double, or to a large extent reabsorb the correc- 
tion due to final state interactions. 

It is clear that in order to obtain a calculation 
of K ^ TTTT with this method one needs to de- 
termine these two constants - and the necessary 
information is not yet available. On the other 
hand, it is important to emphasize the advan- 
tages of this method with respect to CHPT, or 
to the Omnes improved version of it: here one 
does not need to introduce a large number of con- 
stants, but only needs two inputs. The relation 
between these constants and the K ^ tttt ampli- 
tude is exact and it is provided by the solution 
of the dispersion relation. The accuracy of the fi- 
nal result will depend on the accuracy with which 
the two constants are known and the accuracy in 
the solution of the dispersion relation. As far as 
the latter is concerned the prospects of reducing 
it are quite good: the tttt phase shifts, which are 
dominating ]To| , are known with a remarkable ac- 
curacy |l6) . The numerical analysis performed so 
far will be improved by taking into account 
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the effects of the inelastic contributions in the s 
channel [p3| . 

4. CHPT RESULTS ON QUARK MASS 
DEPENDENCE 

In CHPT the dependence on the quark masses 
is explicit: the low energy constants of the ef- 
fective lagrangian are quark- mass independent. 
Therefore, once these are known, one can vary 
the quark masses and see how any given quan- 
tity depends on them. This exercise is particu- 
larly interesting from a lattice perspective, since 
the physical up and down quark masses are still 
far from the reach of today's calculations, and 
one may wish to learn beforehand what to ex- 
pect, when those values are approached. To make 
this exercise worthwhile it is important that one 
goes beyond next-to-leading order: only in this 
manner can one check how fast the series is con- 
verging and where one should stop believing the 
result of the chiral expansion. In the last few 
years there have been a few detailed numerical 
analyses at next-to-next-to-leading order both in 
the S'f/(2) X SU{2) and in the SU{3) x SU{3) 
framework, and I will briefly review the results 
concerning the quark mass dependence. 

4.1. Chiral SU{2) x SU{2) 

In Ref. the low energy constants appearing 
in the tttt scattering amplitude were determined 
by comparing the chiral representation thereof to 
a phenomenological one that follows from Roy 
equations. The numerical solution of these 
provided the necessary input for the determina- 
tion. It has to be stressed that in this procedure, 
only the constants that govern the momentum de- 
pendence of the amplitude could be determined 
- those that give the quark mass dependence had 
to be fixed with theoretical estimates, and as such 
are subject to larger uncertainties. In Fig. || the 
dependence of the two S wave scattering lengths 
on the pion mass is shown. The scattering lengths 
are divided by their leading order value. The fig- 
ure shows that large deviations from leading order 
occur rather early, particularly for the / = case. 
This is an effect of the particularly strong attrac- 
tive interaction in the 7 = channel, as opposed 
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Figure 3. Dependence of the two S wave scatter- 
ing lengths on the pion mass. The ratios Ri^s are 
the scattering lengths normalized to the leading 
order value. 



to the much weaker and repulsive one in the 1 = 2 
channel. The plot shows that lattice calculations 
of these quantities should be made down to rather 
small values of the quark masses if a reliable ex- 
trapolation to the physical value has to be made. 
In comparison, quantities that do not involve two 
pions in the final state, like the pion mass and 
decay constant, show a much milder dependence, 
see Fig. ^. There are calculations of the pion scat- 
tering lengths already available in the literature, 
particularly for the 1 = 2 channel which is consid- 
erably easier. Unfortunately they are all done at 
rather high quark masses. The two most recent 
ones are in Ref. |2^], where reference to earlier 
calculations can also be found. 

The strong quark mass dependence of the 7 = 
scattering length sounds as a warning also for fu- 
ture direct calculations of the K ^ tttt amplitude. 
In that case the FSI are even stronger than for the 
TTTT amplitude at threshold (at s = Mj. the phase 
in the 7 = channel is around 40°), and the quark 
mass dependence should presumably be similarly 
strong ~ implying again the need for simulations 
at rather small values of quark masses, before a 
reliable extrapolation to the physical value can be 
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Figure 5. Dependence on the strange quark mass 
of the meson masses, from ||l7| 

made. 



4.2. Chiral SU{3) x SU{3) 

The calculations of quantities at next-to-next- 
to-leading order in SU{3) x SU{3) CHPT are 
much more complicated, and the convergence 
of the series is certainly slower than for chiral 
SU{2) X SU{2). Recently, an analysis of Kf,4 de- 
cays at this order has been performed In 
this work all the low energy constants of order 
have been estimated with resonance satura- 
tion, whereas the O(p^) constants have been fit- 
ted by comparing to the data on Kei form factors. 
Again, having fixed all the constants, the authors 
were also able to consider the quark mass depen- 
dence of various quantities, and particularly of 
meson masses and decay constants. These are 
shown in Fig. || and |^. The results of this paper 
are puzzling as far as the convergence of the se- 
ries for masses and decay constants are concerned: 
the authors find that the 0{p^) contributions are 
systematically larger than the 0{p'^) ones, as one 
can see from Figs. ^ and ^ It is also surprising 
how strongly the decay constants depend on the 
strange quark mass, including the one of the pi- 
ous. The reader is referred to the original paper 
1l7| for a discussion of the uncertainties involved 



in the determination of the various low energy 
constants - as stressed above, those governing the 
quark mass dependence of meson masses are the 
most difficult ones to control. In view of this, 
these plots should rather be viewed as an indica- 
tion of the qualitative behaviour, rather than pre- 
cise quantitative statements. It would of course 
be very interesting if lattice calculations could 
contribute and calculate the dependence on the 
strange quark mass of various meson masses and 
decay constants. 

5. SUMMARY 

In this talk I have reviewed and critically dis- 
cussed the recent literature on final state interac- 
tions in K — !• TTTT, and a few examples of quark 
mass dependence of various quantities as calcu- 
lated in CHPT at next-to-next-to-leading order. 
I believe that the few issues discussed here are 
good examples of how important the exchange 
of information between the lattice and the chiral 
communities is at present - and of how important 
it will become in the future. 
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Figure 6. Dependence on the strange quark mass 
of the meson decay constants, from [0 
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